Measurement of Two-Qubit States by a Two-Island Single Electron Transistor 
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We solve the master equations of two charged qubits measured by a single-electron transistor 
(SET) consisted of two islands. We show that in the sequential tunneling regime the SET current 
can be used for reading out results of quantum calculations and providing evidences of two-qubit 
entanglement, especially when the interaction between the two qubits is weak. 



Quantum information processing in solid state nanos- 
tructures has attracted wide spread attention because of 
the potential scalability of such devices. Within this con- 
text, quantum measurement in mesoscopic systems is a 
crucial issue and is being carefully analyzed both experi- 
mentally US Hi and theoretically |I B S H H Efl IH 
Ir3 | , so that proper measurements can be designed to ex- 
tract the maximal amount of information contained in a 
solid state qubit (or qubits). One prominent example is a 
single-electron transistor (SET), whose current is partic- 
ularly sensitive to the charge degrees of freedom through 
gate potential variations on its central island(s). Indeed, 
with a radio-frequency SET, electrons can be counted at 
frequencies up to 100 MHz £|, so that if the states of a 
qubit can be distinguished by charge locations, an SET 
can be used to measure the qubit states. 

Recently, two-qubit coherent evolution and possibly 
entanglement have been observed in capacitively coupled 
Cooper pair boxes [l^. The realization and detection 
of two-qubit entanglement are crucial milestones for the 
study of solid state quantum computing. In this Letter 
we study a novel scheme for the quantum measurement of 
two charge qubits (N = 2) , which can be extended to the 
detection of moderately larger number of qubits (N>2). 
Specifically, the target qubits being constantly measured 
are double dot charge qubits whose states are the 
different spatial distributions of the excess electron on 
the double dot. The quantum detector is a two- island 
SET (N = 2), with each island coupled to a qubit ca- 
pacitively, as illustrated in Fig. ffl Our objective is to 
demonstrate the capability of this two-island SET in de- 
tecting and differentiating two-qubit quantum states. In 
particular, we develop a master equation formalism from 
microscipic Hamiltonian to describe the readout current 
of the SET in its sequential tunneling regime. Under 
the condition that the relaxation time of SET current is 
sufficiently long compared to the period of qubit oscil- 
lations, we clarify three major issues regarding the ca- 
pability of the two-island SET layout: whether the two- 
qubit eigenstates {|00), |01), |10), and|ll)} can be distin- 
guished; whether entangled states and product states can 
be distinguished; and whether Zeno effect can be seen in 
the two qubits. 

The Hamiltonian for the combined two qubits and the 
two-island SET can be written as follows: 
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FIG. 1: Qubits are capacitively coupled to a two-island SET, 
which acts as a charge detector. N(>2) qubits are arranged 
between source and drain. 
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FIG. 2: Electronic states in the detector. 



where -ffqb, H sot , and £fi n t are the Hamiltonians of the 
two qubits, the SET, and the interaction between the 
qubits and the SET, respectively. Hqb describes the 
two interacting (left and right, as illustrated in Fig. 1) 
qubits, each consisted of two tunnel- coupled quantum 



dots (QDs) and containing one excess charge 
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where SIl(£Ir) and AL(t)[A R (t)] are the inter-QD (but 
intra-qubit) tunnel coupling and energy difference in the 
left (right) qubit. Here we use the spin notation such that 
<J ax = a^ba + 6+ a a and a az = a^a a ~ tf a b a (a = L, R), 
where a a and b a are the annihilation operators of an 
electron in the upper and lower QDs of each qubit. J 
is a coupling constant between the two qubits, originat- 
ing from capacitive couplings in the QD system 1 f) 
and 1 |) refer to the two single-qubit states in which the 
excess charge is localized in the upper and lower dot, 
respectively. A a (a = L,R) are bias gate voltages ap- 
plied on the qubits, which can be used to tune the qubit 
energy splittings and are used for the manipulation of 
these charge qubits during quantum calculations (lT| . 
The SET part of the Hamiltonian H sct is written as: 
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Here Ci LS {ci RS ) is the annihilation operator of an electron 
in iith (jflth) level (il(*,r) = in the left(right) 

electrode, rfx, s (rfR s ) is the electron annihilation operator 
of the left (right) SET island, s e {T> 4} is the electron 
spin, and n as = <it s <i as is the number of electron on each 
island. Here we assume only one energy level on each 
island. Vl s (Vr s ) and Vms are the tunneling strength 
of electrons between left (right) electrode and the left 
(right) island and that between the two islands. Ul(U r ) 
is the on-site Coulomb energy of double occupancy in the 
left (right) island. Finally, the interaction between the 
qubits and the SET, described by i?; n t, are capacitive 
couplings between the qubits and the two SET islands: 
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Consequently, the energy level of an SET island is raised 
by eC? nt Ct/C%/{ClC? nt + C%[Ct + C? nt ]) if the 

charge in the corresponding qubit is located in the lower 
QD (3 • The electronic states of the qubits also influence 
the tunneling rates T a (E) = 2i:p a (E)\V a {E)\ 2 (p a are 
densities of states of the electrodes). If we define {\A) = 
| ||), \B) = | |T>, |C) = | Tl), \D) = | TT>}, the tunneling 
rates have the relations; T\ — T B < T^, = and T% = 

Now we can construct the equations of the qubits- 
SET density matrix elements governed by the above- 
mentioned Hamiltonian at T = 0, following the pro- 
cedure developed by Gurvitz The possible elec- 
tronic states in the detector are shown in Fig. [2J The 
method is applicable as long as the energy-levels of the 
islands is inside the chemical potential p^ of the left elec- 
trode and fiR of the right electrode, and the tunneling 



rates are much smaller than the difference M£ — /Ur, i.e. 
Ml - PR > {r L , V M } 0. We consider the following 
two transport processes separately. The first case is when 
the double-occupied states are inside the range of pl and 
p R and all electronic states in Fig. 2 take part in the tun- 
neling (finite U model). The second case is when double 
occupancy of electrons [(e)- (i)] is prohibited (infinite U 
model) . Experimentally, these two cases are interchange- 
able by tuning applied island gate voltages |l5j . 

The wave function of the qubits-SET system can 

be expanded over the electronic states of the qubits and 
the island states of the SET shown in Fig. [21 Assum- 
ing that there is no magnetic field and the tunneling is 
independent of spin, after a lengthy calculation, we ob- 
tain 352 equations for density matrix elements Pu\ Z u 2 {t) 
(ui , U2 indicate quantum states of the detector (Fig. 121 
and Zi,Z2 = A, B, C, D are those of the qubits) as |lq |: 

-2^p^-i^ R {p^-pi B )-in L {p c J-p^) 
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(fP = 2^EZ-EZ + E? L+ E? n -J c + J D }-T«')pC D 

+r L '(p^ T +P^), (5) 

where J A = A L + A R + J, J B = A L - A R - J, J c = 
-A l +Ar-J, J d = -Al-Ar+J, Ef L = Eg = E dL +E± t , 

E d L = E d L = E d L - E Lv E d R = E d R = E d a + ER lt ' 

E d a = E d„ = E dR~ E *f r "' = in infinite U model 
and r Q = r Q in finite U model. The readout current 
I(t) = eN R (t) can then be written as [f| 

7Y-A \ ^ f"pi?r zz _l zz _i_ zz i n zz i zz | zz 

— 2-*i ^ CC T ' "eel ' Pd^d^ Pd^d^ Pd n d^ Pd u d u 

z=A,B,C,D 

+ 2 P f f +p^+p^ l +2(p^+pil l )}+2T R 'pli}. (6) 



For simplicity we consider two identical qubits, with 
Ed L =E dR and E int = E(^.=E^ t . 

We monitor the onset of the readout current to ex- 
tract information of the qubit states. The current be- 
gins to flow at t = and after a transient region sat- 
urates to a steady state value. In the meantime, the 
qubits oscillate with frequencies + A^ . The inter- 

action with the dissipative current degrades the coherent 
oscillations and makes the charge distribution uniform 
in the qubits at t — > oo. Conversely, in the absence of 
the qubits, the current saturates around t ~ where 
Y = T L T R /(T L +T R ), while the qubit charge oscillations 
modify the SET current through an effective gate poten- 
tial on the islands. Figure shows the time-dependent 
current characteristics of the infinite U model near t ~ 0. 
At small t state \A) suppresses the current the most while 
state \D) the least. The measurement time t m that is 
required to resolve the states of qubits is estimated as 
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FIG. 3: Time dependent readout current characteristics of 
the infinite U model for \A) = | U), \B) = I IT), \C) = | U), 
|D) = | TT) as initial states (t = 0), where Q L = = 0.75r, 

14/ = o.sr, £ ; L nt = b» = car, = = r« = rg = o.sr, 
ic=r£ = r|=rg = i.2r. (a) j = o.ir, (b) j = r. 
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FIG. 4: Time dependent readout current characteristics 
starting from :(a) \B), (b) singlet state in the infinite U model 
for weak measurement case (-Ejnt = 0.2r < Vm = 0.5T) as a 
function of V s = Vg = V^. Parameters are the same as those 
in FigEl 



f- 1 ~min{£' illt ,r^-r^} (~ 0.5 -1 r in Fig. EJ. The rela- 
tive magnitude of the current changes after the coherent 
motions of qubits (t > l/f2 a ). Thus the SET current 
can be used to distinguish the four product states during 
t m < t < 1/Oq.. If the coherent oscillation of the qubits 
remains after t > T -1 , as in the present model 0], we 
can discuss the quantum states of qubits using the steady 
current formula (t — > oo) through the SET without the 
qubits 0: 
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where = E^ L — E^ R is the energy difference of the 
two islands. If Vm 3> I\ e<j, Q a > the coupling between the 
two islands is strong and the current mainly reflects the 
bonding-antibonding state in the detector, which is not 
suitable for qubit measurements. We thus focus on the 
regime of V M <Q a ,T. Since -Ej R = Eg- Eg = 
and Eg-E? = E?-Eg = 2E- mt , the different ef- 
fects between \A) and }D) and that between \B) and \C) 
come from the differences in the tunneling rates. More- 
over, the difference of \A) and \D) from \B) and \C) be- 
comes obvious in the E^ > Vm region. Thus we call 
Eint > Vm strong measurement regime, where the four 
product states can be distinguished, in contrast to the 
weak measurement regime of E mt < Vm ■ 

We can distinguish the current of pure entangled states 
and that of pure product states by changing bias volt- 
ages Vg — A Q in the regime of J/T <C 1, where the cur- 
rent depends on the change of qubit oscillation frequency 
(~ + ). Figure a) shows the current corre- 

sponding to the qubit \B) state in the weak measurement 
regime of the infinite U model. We also obtained similar 
results for the other product states \A), \C), and \D). In 
contrast, the readout current for a two-qubit entangled 
state is more uniform compared with the product states 
as entangled states generally have less distinct charge dis- 
tributions. For example, the density matrix elements for 
a singlet state (|U) - |IT))/V2 = (|C) - \B))/y/2 of two 



free qubits (H- mt = 0) satisfy p UB + p CG - p BC - p CB = 
(Al = Ar), which suggests that entangled states such 
as the singlet state are less effective in influencing the 
readout current. We believe this ineffectiveness is re- 
lated to the fact that logical states encoded in entangled 
states are less susceptible to environmental decoherence 
(l8j . Indeed, the readout current of this entangled state 
is found to be uniform as shown in Fig. |l£b). We ob- 
tained similar results for the other Bell states, and there 
is no significant difference between the infinite U model 
and the finite U model in the weak measurement regime. 

In the strong measurement regime (E[ nt > Vm), the 
current is more sensitive to the charge distributions in 
the qubits, and there are differences between the infi- 
nite U model and finite U model. We can distinguish 
the four products more easily through the SET current, 
as shown in Fig. Ela)-(d). However, currents for the 
entangled states in the infinite U model show several 
similar peaks that reflect the qubit oscillations and can- 
not be easily distinguished from the product states. On 
the other hand, the finite U model shows distinct uni- 
form structure compared with the current of the product 
states [Fig.[!Je) and (f)]. This shows that, in the finite U 
model, redistribution of the electrons through the two is- 
lands of the detector is energetically favorable under the 
rather uniform electric field generated by the entangled 
qubits. Figure EJa) shows that the concurrence (a mea- 
sure of entanglement [l9l | derived from reduced density 
matrix of two qubits after tracing over the detector com- 
ponents) of the two qubits disappears quickly in the cases 
of strong measurement. While the coherence quickly de- 
grades, we can see the emergence of the Zeno effect, in 
which a continuous measurement slows down transitions 
between quantum states due to the collapse of the wave- 
functions into observed states 0,^3- F° r instance, Fig. El 
(b) shows that, as E- mt increases, the oscillations of den- 
sity matrix elements of the qubits (e.g. p DD ) are delayed, 
which is a clear evidence of the slowdown described by 
the Zeno effect in the two qubits. 

Our numerical results above are applicable to a wide 
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FIG. 5: Time dependent readout current characteristics in 
the finite U model (U = 2r)for strong measurement case 
(E int = 0.8r > Vm = 0.5r) as a function of V g = Vg = V* . The 
initial states are (a) \A), (b) \B),(c) |C),(d) \D), (e) triplet 
state, (f) singlet state. Parameters other than Ei B t are the 
same as those in Figl3*l 



coherent oscillations of the qubits when the gate bias 
changes between V g L = V g R and V g L = — V g R . 

Since the detection scheme discussed here is based 
on measuring small current differences in the transient 
regime, it is important to analyze whether the present 
day technology can achieve the necessary sensitivity. The 
state of the art technology allows the measurement of 1 
pA current with dynamics in the GHz frequency range 
with repeated measurement techniques 0, 0, |2(| . Ac- 
cording to our Figs. 3-5, our scheme requires measuring 
a 0.1 pA current that changes in the nanosecond time 
scale (assuming a T in the order of 100 MHz, a reasonable 
figure because E- mt would be in the order of 100 MHz if 
all capacitances are 100 aF), which is at the edge of the 
current measurement technology. Thus, with a similar 
design of repeated measurement 0, 0, |2(J , our detection 
lio 1 
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FIG. 6: (a) The concurrence of the singlet state, (b) Example 
of Zeno effect: oscillation of p DD (t) is delayed, where the 
initial state is \D) state (p DD (0) = 1). Similar effects can be 
seen in other initial states. Parameters are the same as those 
in FigEl 



range of pure product and entangled states. For exam- 
ple, in the entangled states cos0| Ti) + e iV sin#| ||), we 
found that the uniformity of the readout current holds 
approximately up to \9 ± 7r/4|, \<p\ < it/12. The pure 
entangled states are more robust beyond the spatial dis- 
tribution of the wave functions. Although the product 
statesnorf^Icosf^e-'^l t) Q +sin(^)e 1 ^ | |) Q ] seem 
to have similarly uniform wave functions when 9l = ±9r 
and (fL — ±V?-R — 0j n (compared to the entangled states 
mentioned above) , the corresponding currents reflect the 



scheme should be experimentally feasible in the near fu- 
ture. 

In conclusion, we have solved master equations and de- 
scribed various time-dependent measurement processes 
of two charge qubits. The current through the two- island 
SET is shown to be an effective means to measure results 
of quantum calculations and entangled states. 

We acknowledge discussions with N. Fukushima, S. Fu- 
jita, and M. Ueda. 
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